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Abstract 

The existence and construction of common invariant cones for families of real 
matrices is considered. The complete results are obtained for 2x2 matrices 
(with no additional restrictions) and for families of simultaneously diagonal- 
izable matrices of any size. Families of matrices with a shared dominant 
eigenvector are considered under some additional conditions. 
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1. Introduction 

The theory of nonnegative matrices, and more generally of matrices that 
leave invariant a convex, closed, pointed, solid cone, is classical; we mention 
here the books [1, 2] among many others; see also [3] for a review of many 
results, including recent ones, and extensive bibliography. More generally, 
real matrices that leave invariant a convex, closed, pointed, solid cone, have 
been studied in [4, 5, 6, 7, 8, 9]. A complete characterization of such matrices 
in terms of spectral structure was obtained in [5] . An interesting application 
to the multiple agents randezvous problem is given in [10]. 
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Recently, several works appeared studying matrices having common in- 
variant convex, closed, pointed, solid cones. These works have been mo- 
tivated primarily by applications in Glass networks [11] and joint spectral 
radius [12, Theorem 1]. Glass networks are continuous-time switching net- 
works used to model gene regulatory networks and neural networks; see [11] 
and references there for an in depth discussion on Glass networks. 

The paper [11] actually served as a motivation for the current paper. We 
develop here results on matrices having common invariant cones. The auxil- 
iary Section 2 contains necessary notions and definitions, in particular that of 
a proper cone and a dominant eigenvector. In Section 4, a full description is 
given of families of 2 x 2 real matrices having common invariant proper cones. 
As it turns out even in this case the characterizations are rather involved, 
and the proofs not immediate. Some partial results (for pairs of diagonaliz- 
able but not simultaneously reducible matrices) in this venue were obtained 
in [11]. Our approach is based on the description of all invariant cones for a 
single 2x2 matrix given in Section 3. In spite of its elementary nature, we 
did not find this description in the literature, and include it for the sake of 
self containment. Section 5 contains the existence criterion for (and actually 
a construction of) a common invariant cone of a family of simultaneously 
diagonalizable matrices, while Section 6 provides some sufficient conditions 
for such a cone to exist when the matrices share the dominant eigenvector. 
Finally, Section 7 consists of several examples illustrating both the results 
obtained and their limitations. 

2. Preliminaries and definitions 

Let M be the field of real numbers, M"^ the set of real n-component column 
vectors, and R™^'^ the set of real mxn matrices. All matrices in the present 
paper are assumed to be real, unless explicitly stated otherwise. A set /C C M" 
is a cone if a/C C /C for all scalar multiples a > 0. A cone /C is said to be 
proper if /C + /C C /C (so that /C is convex), closed, pointed {K fl — /C = {0}) 
and solid (the interior of /C is nonempty). 

For X being a subset of M" or we denote by ConeX the smallest 

convex cone containing X and say that X generates ConeX. Of course. 
Cone X is nothing but the set of all (finite) linear combinations of elements 
of X with non-negative coefficients. A cone having a finite generating set is 
called polyhedral. Polyhedral cones are always closed. 
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For a square matrix A, by the degree of its eigenvalue A in this paper we 
understand its muhiphcity as a root of the minimal polynomial of A (that is, 
the size of the largest block, in the Jordan canonical form of the matrix, cor- 
responding to the eigenvalue A). We will denote the eigenvalues of an n x n 
matrix A by Xi{A), . . . , Xn{A) (or simply by Ai, . . . , A„ if the choice of the 
matrix is clear from the context), always taking p{A) = Ai provided that the 
spectral radius p{A) of A is an eigenvalue. Wc will call the respective eigen- 
vector (eigenspace) the dominant eigenvector (jcsp., dominant eigenspace) of 
A. In case when an eigenspace is one dimensional, we will (naturally) call it 
an eigenline. We will also use the term eigenray for each of the two rays into 
which an eigenline is partitioned by the origin. Finally cr{A) will be used to 
denote the set of all eigenvalues of A. 

A cone IC C M" is said to be invariant under A G M"^" if Ax G /C for every 
X E K,. The following remark is trivial, but will be useful in our analysis. 

Remark 1. A cone /C = Cone {vi, ... ,Vm} is A-invariant if and only if 
Avj G JC for j = 1, 2, . . . , m. 

The following result was proved by Vandergraft [5]. 

Theorem 1. A G IR"^" has an invariant proper cone if and only if 

(i) The spectral radius p{A) G cr{A), and 

(ii) degAi(A) > degXi{A) for every eigenvalue Xi{A) with \Xi{A)\ = Xi{A). 
If conditions (i)-(ii) hold, then also 

(iii) Any A-invariant proper cone contains a dominant eigenvector of A. 

For spectral criteria for existence of polyhedral proper invariant cones see 
[6, 9]. 

We will be using the term Vandergraft matrices for real matrices satisfying 
conditions (i) and (ii) of Theorem 1, denoting the set of all such nxn matrices 
by E^^". 

3. Invariant proper cones for 2x2 matrices 

It is very easy to characterize matrices in M^^^. Namely, condition (i) of 
Theorem 1 is equivalent to 

(traceA)^ > 4detA, traced > 0, (1) 
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the first inequality in (1) meaning simply that the eigenvahics of A are real 
while the second inequality guarantees that the one with the bigger absolute 
value is non-negative. Since condition (ii) then holds automatically, a 2 x 2 
matrix A is Vandergraft if and only if it satisfies (1). 

Conditions (1) hold, in particular, when both eigenvalues Ai,A2 of A 
are non-negative. Description of all 74-invariant proper cones in this case is 
given by the following two theorems, dealing with diagonalizable and non- 
diagonalizable matrices A separately. Of course, in the former situation only 
the case Ai 7^ A2 is of interest, because otherwise A is a scalar matrix which 
leaves every cone invariant. 

Theorem 2. Let a 2x2 matrix A be diagonalizable, with Ai > A2 > 0. Then 
a proper cone fC (ZM.^ is A-invariant if and only if it contains an eigenvector 
of A corresponding to Ai and its interior does not intersect the eigenline of 
A corresponding to A2. 

Proof. "Only if part. An A-invariant proper cone JC must contain an eigen- 
vector of A corresponding to Ai, as follows from Theorem 1, part (iii). Denote 
this vector by Ui and suppose for a moment that there is an eigenvector U2 
of A corresponding to the eigenvalue A2 and lying in the interior of /C. Then 
for sufficiently large M > also —Ui + Mu2 G /C, and for all n = 1, 2, . . ., 



Letting — 00, from the closedness of /C we conclude that —ui e /C. This, 
however, contradicts pointedness of /C. 

"If part. Any proper cone in is generated by two linearly independent 
vectors: A — Cone{^;i, i'2}. The conditions imposed on /C mean that, after 
appropriate scalings, its generating vectors can be written as 



where x > 0. (Here ui, U2 are eigenvectors corresponding to Ai, A2, respec- 
tively.) Then 



(ArM)"(-Mi + MM2) 



ui + M(A2/Ai)"m2 e IC. 



Vi= Ui^- U2, V2 = Ui— XU2, 



Avi{^ AiWi \2U2) 



xXi + A2 



V2 e K, 



l + x 



1 + x 



1 + a; 

A-invariance of /C therefore follows from Remark 1. 




Ai x\2 



l+x 



V2 e /C. 



□ 
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Let now A e R^^^ be non-diagonalizable. Then, for any v e R^, 

Av ^ Xv + xu, (2) 

where A is the eigenvalue of A, u is its (arbitrarily fixed) eigenvector, and 
X e R. We will say that v is positively /negatively associated with u {relative 
to A, if there is a need to mention the matrix explicitly) if in (2) ±x > 0. 
Observe that x = if and only if v belongs to the eigenline of A, that is, is 
a scalar multiple of u. 

Of course, v is positively associated with u if and only if —v is negatively 
associated with u if and only if —v is positively associated with —u. Geo- 
metrically speaking, the plane is partitioned by the eigenline of A into 
two open half-planes; one consisting of vectors positively associated with u, 
and the other of vectors negatively associated with u. 

Theorem 3. Let A G R^^^ be a non-diagonalizable matrix with the eigen- 
value A > 0. Then a proper cone JC is A-invariant if and only if it is given by 
K, — Cone{M, v}, where u is an eigenvector of A and v is positively associated 
with u relative to A. 

Proof. "IF part. Since A > 0, from (2) it follows that Av e Cone{u,v}, 
because a; > 0. Obviously, Au = Xu also lies in Cone{u,v}. The desired 
result now follows from Remark 1. 

"Only if part. Let a proper cone /C be A-invariant. Due to Theo- 
rem l(iii), there is an eigenvector of A lying in /C. Denoting it by u, observe 
that vectors negatively associated with u cannot lie in fC. Indeed, if A = 
and (2) holds with x < 0, then 

V & K. =^ —u e /C, 

which contradicts the pointedness of JC. For A > 0, (2) implies 

A'^v = X'^v + nxX^'-^u, n = l,2,.... 
Consequently, if v e /C and x < 0, then 

-u = lim ^A^-"A"t; e /C 

n^oo n \x\ 

— once again, a contradiction with the pointedness of /C. 

Since in every neighborhood of u there are vectors negatively associated 
with it, u cannot lie in the interior of /C. Thus, it must be one of its generating 
vectors. The other generating vector v, being linearly independent with u, 
must be positively associated with it. So, /C indeed is of the desired form. □ 
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Corollary 1. For non-diagonalizable Vandergraft 2x2 matrices, the domi- 
nant eigenvector lies on the boundary of their invariant proper cones. 



As follows from Theorem 2, for diagonalizable 2x2 matrices with posi- 
tive eigenvalues the dominant vector can lie both in the interior and on the 
boundary of their invariant cones. 

We turn now to the remaining case of matrices A with negative determi- 
nants. Denote the eigenvalues of A by Ai(> 0) and A2(< 0), and let ui,U2 
stand for the respective eigenvectors. 

Theorem 4. Let A e M^''^ and det A < 0. Then a proper cone IC G M. is 
A-invariant if and only if it can he represented as K. — Cone{i;i, ^2}; where 

Vj ^ui + CjU2 (j = 1, 2) (3) 

and 

ci>0, C2<0, ^<^<^ (4) 

A2 C2 Ai 

Proof. An A-invariant pointed cone cannot contain eigenvectors of A corre- 
sponding to a negative eigenvalue. Thus, all vectors f G /C (in particular, its 
generators) in their expansion along the eigenbasis {mi,M2} have the same 
sign coefficients corresponding to ui. Switching from ui to —ui if needed, we 
may without loss of generality suppose that these coefficients are positive. 
Scaling Vi and V2 if necessary, we arrive at (3). Yet another change (from U2 
to —U2, or flipping Vi with V2) allows us without loss of generality suppose 
that Ci > C2. 

On the other hand, for Vj given by (3) we have 

AVj = Xi{ui + ^CjU2). 

Ai 

Consequently, Avj lie in the cone /C if and only if the numbers X2Xi^Cj lie in 
[c2,Ci]. This is equivalent to (4). □ 

Corollary 2. Let A be a 2x2 Vandergraft matrix with negative determinant. 
Then any dominant eigenvector of A lies in the interior of those A-invariant 
proper cones that contain the eigenvector. 

Note that conditions (4) are consistent if and only if det A < and 
trace ^4 > 0, which of course agrees with (1). If this is indeed the case, for 
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every non-zero vector v different from the eigenvectors of A there exist A- 
invariant proper cones /C with v being one of the generators. The second 
generators of these cones form yet another convex cone, described by (4) 
with one of cj being determined by v and the other serving as a parameter. 
The latter cone degenerates into a single ray if and only if trace A = 
(equivalently: A^ is a scalar multiple of the identity), when necessarily ci = 

-C2. 

It is very easy to produce directly an A-invariant cone with arbitrarily 
chosen generator v for any 2x2 matrices A with 

det A<0, trace A>0. (5) 

Lemma 5. Let A E M^^^ satisfy (5). Then K,y := Cone{'t;, ^It;} is A- 

invariant for any v E'M? , v ^ 

Of course, Ky is proper if and only if v is not an eigenvector of A. 

Proof. Indeed, is generated by v and Av. The first of these generators is 
mapped by A into ICy by construction, and 

A{Av) = A^v = (trace + (-det A)^; e /C^ 

due to the Cayley-Hamilton theorem. □ 

This simple observation will become useful in the next section. 

4. Common invariant cones for families of 2 X 2 matrices 

Let A — {Ai, . . . , An] be a finite family of 2 x 2 real matrices. An A- 
invariant proper cone by definition is A^-invariant for all j = 1, ... , n, and in 
order for that to be possible each of the A^ 's must be a Vandergraft matrix. 

In particular, matrices of the form cl with c < preclude the existence 
of ^-invariant proper cones. On the other hand, presence (or absence) of 
matrices cl with c > in ^ is irrelevant. All such matrices (if any) can be 
deleted from A but may as well be left intact. 

We first consider the case when all the matrices Aj share a dominant 
eigenvector u. If several of them are non-diagonalizable, we will say that they 
have the same orientation if the sets of vectors positively associated with u 
relative to these matrices coincide (of course, the sets of vectors negatively 
associated with u then coincide as well). This happens if and only if in a 
basis containing u the off diagonal elements of these matrices are all of the 
same sign. 
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Theorem 6. Let A — {Ai, . . . , A^} be a family of 2x2 Vandergraft matrices 
sharing the same dominant eigenvector u. Then there exists an A-invariant 
proper cone K, if and only if either 

(i) all Aj are diagonalizable, and those of them (if any) which have 
det Aj < 0, trace Aj — are scalar multiples of each other, or 

(ii) all Aj have non-negative determinants, and those of them which are 
not diagonalizable (if any) have the same orientation. 

Proof. "If part, (i) If all Aj are diagonalizable and have non-negative deter- 
minants, the result follows from Theorem 2: any proper cone /C containing u 
and sufficiently narrow to avoid all the eigenvectors of Aj corresponding to 
their second eigenvalue will do the job. 

Suppose now that some of Aj have negative determinants; relabel them 
by Ai, . . . , Ak. Consider /C = Cone{u, v, Aiv, . . . , Akv}, where v is a vec- 
tor different from u but so close to it that K. does not contain the non- 
dominant eigenvectors of A^+i, . . . ,An and AiAj {i,j = 1, . . . , k). The lat- 
ter products are all Vandergraft matrices with positive determinants and, 
under conditions imposed, also diagonalizable (here, the hypothesis that 
all matrices Aj with negative determinants and zero traces are multiples 
of each other, is crucial). Hence, /C is invariant under Ak+i,...,An and 
AiAj (i,j — 1, . . . , Ac) as in the previous part of the proof (see Theorem 2). 
In particular, AiAiv, . . . , AiA^v G /C for all i = 1, . . . , fc. Since AiV G /C 
{i = 1,2, ... ,k) by the construction of /C, in fact all the generators of IC 
are mapped by Ai, . . . ,Ak into /C, so that JC is invariant also under Ai for 
i — 1, . . . ,k. 

(ii) There is no need to consider the case when all Aj are diagonalizable, 
because it is covered by (i). Supposing that non-diagonalizable matrices 
are present in A, relabel them by ^41, . . . ,Ak. Choose a vector v positively 
associated with u relative to Ai, under the conditions imposed it will be 
positively associated with u also relative to A2, . . . , Ak. By Theorem 3, K, — 
Cone{u,v} is Aj-invariant for j = l,...,k. Moving v sufficiently close to 
u in order to avoid the non-dominant eigenvectors of A^+i, . . . , A^, we will 
make JC invariant with respect to all Ai, . . . , An. 

"Only if part. In cases different from (i)-(ii) the family A contains either 
(iii) two linearly independent matrices with negative determinants and zero 
traces, or (iv) two non-diagonalizable matrices with different orientation, or 
(v) a non-diagonalizable matrix and a matrix with negative determinant. 
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Denote the matrices involved in each case by Ai and Then in case 
(iii) A1A2 is a non-diagonahzable Vandergraft matrix, so that (iii) reduces to 
(v). In case (iv), due to the description given by Theorem 3 the intersection 
of any Ai-invariant proper cone with an y42-invariant proper cone is a ray 
spanned by u, and therefore not proper. In case (v), the non-existence of 
common invariant proper cones follows from the comparison of Corollaries 1 
and 2. □ 

Corollary 3. In the setting of Theorem 6, an A-invariant proper cone exists 
if and only if any two matrices in the family A share an invariant proper 
cone. 

Proof. Indeed, from the consideration of cases (iii)-(v) in the proof of The- 
orem 6 it follows that there exists a pair of matrices in A with no common 
invariant proper cone, whenever conditions (i) or (ii) do not hold. □ 

We now move to the situation when A contains matrices with different 
dominant eigenlines. As it happens, the crucial role is then played by an 
extended family Ai which contains A and all pairwise products (different 
from scalar multiples of the identity) of the matrices in A having negative 
determinants: 

Ai^ AU {AiAj : Ai, Aj e A, det A^, det Aj < and AiAj ^ c/}. 

Of course, A\ coincides with A if the latter consists only of matrices with 
non-negative determinants. 

Theorem 7. Let A be a finite family in R^^^. For an A-invariant proper 
cone to exist it is necessary that 

(i) all elements of Ai are Vandergraft matrices, 

(ii) there are at most two dominant eigenlines corresponding to 
non-diagonalizahle matrices in Ai, and all of them (if there is more 
than one) corresponding to the same dominant eigenline also have the 
same orientation, 

and 

(iii) the dominant eigenlines of matrices in Ai are separated from the non- 
dominant ones. 
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The separation condition (iii) means simply the existence of vectors Vi, V2 
such that the interior of Cone{?;i, ?;2} is free of the non-dominant eigenvectors 
of matrices in Ai while the interior of Conejfi, —V2} is free of the dominant 
eigenvectors of non-scalar matrices. The vectors Vj themselves are allowed to 
be both dominant and non-dominant, but only if as the latter they correspond 
to matrices in Ai with non-negative determinants. 

Proof. An ^-invariant cone fC also is ^i-invariant. This immediately implies 
the necessity of condition (i) . 

According to Corollary 1, an eigenline of a non-diagonalizable 2x2 Van- 
dergraft matrix must be passing through the boundary of any of its invariant 
proper cones. Thus, at most two such eigenlines are admissible. 

If two non-diagonahzable matrices share the eigenhne but have different 
orientation, the intersection of (any pair of) the respective invariant cones 
is an eigenray, due to Theorem 3, and therefore is not proper. These two 
observations settle the necessity of part (ii). 

Finally, if JC is an A- (and therefore v4i) -invariant proper cone, then all 
dominant eigenlines lie in /C U (— /C) while non-dominant eigenlines belong 
to the closure of the complement (to the complement itself, if the respective 
matrix has negative determinant). Thus, (iii) holds. □ 

Suppose now that necessary conditions stated in Theorem 7 hold. Denote 
by U — {ui, . . . ,un} the set of all distinct dominant unit eigenvectors of 
matrices in Ai the directions of which are chosen in such a way that Cone U 
is proper and its interior is free of non-dominant eigenlines (this is possible 
due to (iii)). If there are no such eigenlines (that is, all matrices in A are 
non-diagonalizable), impose instead the condition that Uj for j = 2, . . . ,N 
are positively associated with ui relative to the matrix Ai for which ui is an 
eigenvector (this is possible due to (ii)). This choice is unique up to changing 
the sign of all Uj simultaneously. Relabel also the elements of A in such a 
way that det is negative for i = 1, . . . ,k and non-negative otherwise (with 
the convention that k — if det > for alH = 1, . . . , n). 

For further consideration it is convenient to distinguish between the cases 
when there is none, one, or two dominant eigenlines corresponding to non- 
diagonalizable matrices in Ai- 

Theorem 8. Let A — {Ai, . . . , An} C R^^^ be such that all the elements of 
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Ai are diagonalizable matrices. Under the necessary conditions^ (i), (iii) of 
Theorem 7 and using the notation introduced above, let 

K. = Cone{uj, AiUj : i = 1, . . . , k; j = 1, . . . , N}. (6) 

Then there exist A-invariant proper cones if and only if the cone JC is proper, 
its interior is free of the non-dominant eigenvectors of all matrices in Ai, 
and the edges ofK. are not collinear with the eigenvectors of (i = 1, . . . , k). 

Proof. "Only if part. Any ^4- invariant cone also is ^i-invariant, and thus 
must contain either U or —U. Without loss of generality, let it contain U. 
Then, being invariant under all A^, it must also contain /C. The rest follows 
from Theorems 2 and 4, applied to each of the matrices in Ai- 

"If part. For i = k + 1, . . . ,n, the cone /C contains the dominant eigen- 
vector of Ai (since it is one of the u/s) and the interior of /C does not contain 
its non-dominant eigenvectors. By Theorem 2, JC is invariant under A^. 

Since det AiAm > for all i,m = 1,...,A;, the cone /C for the same 
reasons is AiAm-i'nvaxiant. Consequently, A^A^Uj e /C for alH, m = 1, . . . ,k; 
j = 1, . . . ,N. But AiUj lies in /C by construction. So, all the generators of 
/C are mapped into IC hy Ai, . . . , A^. It remains to invoke Remark 1. □ 

Theorem 8 shows that necessary conditions stated in Theorem 7 in gen- 
eral are not sufficient. For k = 0, however, /C coincides with Conef/, and 
the latter is ^-invariant already under the conditions of Theorem 7. The 
situation therefore simplifies as follows. 

Corollary 4. Let A = {Ai, . . . , An} be a family of diagonalizable 2x2 
matrices with non-negative determinants. Then in order for an A-invariant 
proper cone to exist it is necessary and sufficient that 

(i) all elements of A are Vandergraft matrices, and 

(ii) the dominant eigenlines of matrices in A are separated from the non- 
dominant ones. 

We can now observe the following. 



Condition (ii) holds automatically. 
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Theorem 9. Let A — {^i, ■ ■ ■ , A^} he a family of diagonalizable 2x2 Van- 
dergraft matrices with non-negative determinants. If any four of them (three 
— if there is at most one pair of simultaneously diagonalizable matrices in A) 
have a common invariant proper cone, then there also exists an A-invariant 
proper cone. 

Proof. Indeed, if an ^4- invariant proper cone does not exist, then condition 
(iii) of Theorem 7 fails. But then it is possible to find four matrices in 
A (without loss of generality relabel them by Ai, . . . ,A^ such that, when 
traveling around the origin in a counterclockwise direction, one encounters 
consequently the dominant eigenline of Ai, the non-dominant eigenline of A2, 
the dominant eigenline of A3, and finally the non-dominant eigenline of A4. 
Condition (iii) fails for the set {Ai, y42, A3, ^44}, so that these four matrices 
already do not have a common invariant proper cone. Of course, it is not 
excluded that Ai or A^, coincides with A2 or A4, and then we have an even 
smaller subfamily of A with no common invariant proper cone. If Ai and 
A3 are not simultaneously diagonalizable, the non-dominant eigenline of at 
least one of them will be different from the dominant eigenline of the other. 
Consequently, in this case we can always choose A2 or A4 coinciding with 
Ai or A3. A similar reasoning works if a pair A2,A4 is not simultaneously 
diagonalizable. □ 

We now move to the case of one dominant eigenline corresponding to 
non-diagonalizable matrices. 

Theorem 10. Let A = {Ai, . . . , A„} C R^^^ satisfy conditions (i)-(iii) of 
Theorem 7, with exactly one dominant direction (say, corresponding to ui) 
shared by all non-diagonalizable matrices in Ai- Then there exists an A- 
invariant proper cone if and only if the cone (6) satisfies conditions of The- 
orem 8, and in addition its interior consists of vectors positively associated 
with ui. 

Proof. "Only if part. An ^-invariant cone is invariant under all diagonal- 
izable matrices in Ai- Therefore, it must contain the cone (6). On the other 
hand, it is also invariant under all non-diagonalizable matrices in Ai, so that 
by Theorem 3 Ui must lie on its boundary, and the interior of the cone con- 
sists only of vectors positively associated with Ui. The same is therefore true 
for /C. 

"If part. As in Theorem 8, the cone (6) itself does the job. □ 
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The case of two dominant eigenlines corresponding to non-diagonalizable 
matrices in Ai can be treated along the same hnes. However, a more straight- 
forward (and less computationally consuming) approach also is available. 

Suppose that conditions (i), (ii), and (iii) of Theorem 7 hold, and that Ai 
contains two non-diagonalizable matrices (say, Bi and B2) with non-coUinear 
dominant eigenvectors. Relabel the latter as Ui and U2, choosing the direction 
of Ui arbitrarily, and the direction of U2 in such a way that it is positively 
associated with ui relative to Bi. According to Theorem 3, then either Bi 
and B2 have no common invariant proper cones (if Mi is negatively associated 
with U2 relative to B2), or there are exactly two such cones: /C = Cone{Mi, U2} 
and — /C. 

Theorem 11. For a finite family A = {Ai, . . . , An} of Vandergraft matri- 
ces with exactly two dominant eigenlines corresponding to non-diagonalizable 
matrices in A\, the only possible A-invariant proper cones are ±/C introduced 
above. These cones are indeed A-invariant if and only if: 

(a) all non-diagonalizable matrices in A (if any) with a dominant eigen- 
vector Uj have the same orientation as Bj (j — 1,2), 

(b) for all matrices Aj G A, their dominant eigenvectors lie in fCU (— /C) 
while the non-dominant ones lie outside the interior of K, Li (— /C), 

(c) for A = Aj & A with the eigenvalues Xij > 0, \2j < and the domi- 
nant eigenvector ui + ^U2 G /C, the non-dominant eigenvector must be 
collinear with u\ -\- r]U2, where 

\2j i 

Proof. Indeed, conditions (a)-(c) arc necessary and sufficient for /C (or — /C) 
to be invariant under all matrices in A, as follows by applying Theorem 2- 
4. And, as was observed earlier, no other proper cones can possibly be 
^-invariant. □ 

Remark 2. It follows directly from the proof of Theorem 11 that if in its 
setting every three matrices in Ai ( or any five matrices in A) have a common 
invariant proper cone, then there also exists an A-invariant proper cone. 
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5. Simultaneously diagonalizable matrices 

We now move to square matrices of arbitrary size mxm. but suppose that 
all the elements of the family A — {Ai, . . . , An} under consideration can be 
put in a diagonal form by the same similarity transformation S (note that 
S is allowed to be a complex matrix). This S then diagonalizes all matrices 
from A2 = Cone A, and moreover from the closed algebra ^3 generated by A. 
Denote hy q (< m) the maximal number of distinct eigenvalues for matrices 
in A2- If -Bo is one of the matrices on which this number is attained, 

Bo = Sdiag[biIs^,...,bqIs^]S'^, bj e C, bi^bjiii^ (7) 

then 

Aj = 5'diag[Aij/5i, . . . , \jIs^]S''^ for all Aj e A; here Xij e C. (8) 

Indeed, if at least one of the blocks in the middle factor of (8) were different 
from a scalar multiple of the identity, then the matrix Bq + eAj would have 
more than q distinct eigenvalues for sufficiently small e. (Note in passing, 
though this fact is not needed in what follows, that because of (8) q is also 
the maximal number of distinct eigenvalues of the matrices in a larger set 

Assume there exists an ^-invariant proper cone. Then obviously all prod- 
ucts A^^---A^" [rrii e Z+, the set of nonnegative integers) are Vander- 
graft matrices. Due to the diagonalizability, this requirement amounts to 
maxj I A"^'^ ■ ■ ■ A™" | being attained on some i for which A^^ • • • A^" > 0. For 
every n-tuple (mi, . . . , m„) G , introduce the set 

Q(mi,...,m„) = {ioe{l,2,...,g} : max |A7 • • • A^l = A.'^i • • • A,'^;}. 

Although r2(mi, . . . ,m„) need not be a singleton, we note that there is a 
unique index p = p{mi, . . . , m„) G f2(mi, . . . , m„) for which 

max \big\ = bp. 

ioeii{mi,...,m„) 

Indeed, this follows from the Vandergraft property of the matrix 

and the condition bi 7^ bj Hi ^ j. (Note that wc take X° = I for every square 
matrix X regardless if X is singular or not.) We let P — Up(mi, . . . , m„) 
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where the union is taken over all n-tuples (mi, . . . ,m„) G Z". Permuting 
the columns of 5* if necessary, we may suppose without loss of generality that 
this set is P = {1, . . . , k}, where k < q. 

Theorem 12. In the notation (8) and for k as introduced above, A-invariant 
proper cones exist if and only if 

Kj > for all i — 1, ... k and j — 1, . . . ,n. (9) 

Proof. "Only if part. For an arbitrarily fixed io G {1, . . . , A;}, pick an n-tuple 
mi, ... , rUn such that 

Ai:{---A-; = .max|A----Ar|. 

1=1, ...g 

Then 

for any e > 0, and therefore 

a;:} • • • Ai:;; + e + 5b,, >\X-l^...X2-+e + 5k\, i ^ io 

for 6 > small enough. Having fixed e and 6 (> 0), observe that then for 
any j such that Aj^j 7^ 0, 

I iKl ■■■Kn + ^ + ^K)'K3 1 > I (a;^^ • • • Ar + e + h ^ 7^ ^0 

if the positive integer I is large enough. 

In other words, (A^J • • • A^^ + e + Sbi^yXi^j is strictly bigger (by absolute 
value) than other eigenvalues of 

Bi := {AT' • • • + + (^5o)'>l,-. 

But an ^-invariant cone is also P^-invariant whenever e,d > 0. So, 

Choosing two consecutive values of /, we conclude that in fact Aigj > 0. 

"If part. Denote by L_|_ the (real) linear span of the first si + • • ■ + Sk 
columns of S. Note that since the eigenvalues of Bq corresponding to these 
columns of S are real (see (7)), the first si H — • + Sfc columns of S are real as 
well (or more precisely can be made real if necessary, by (complex) scalings); 



15 



thus C M™". Let us represent R"^ as the direct sum of the subspaces 
and Lc spanned respectively by the real columns of S and by the real and 
imaginary parts of non-real (if any) columns of S. By definition of L+, it 
lies in L^. Moreover, can be written as Lr = + L_, where L_ is also 
spanned by columns of S. 

Choose bases F± in L± consisting of columns of S, and a basis Fc in Lc 
consisting of vectors Ui,Vi e M.^ such that 

AjUi = (Re Xij)ui - (Im Xij)vi, AjVi = (Im Xij)ui + (Re Xij)vi. 

Then of course 

(ArMr---^rK = (ini//^M+(Re/i>,, ^ ' 

where rrij G Z+ and /Xj = A"i' • • • A^". 

Denote by / the sum of all elements in F+, and let /Co stand for the 
smallest ^-invariant convex cone containing F+, / + F_ and / + Fc- The 
span of /Co contains the basis F = F^U F_ U F^ of the whole space M™, so 
that it coincides with M*". In other words, JCq is a reproducing convex cone, 
and therefore it is solid. 

The closure /C of /Co also is a convex solid cone invariant under A. It 
remains only to show that /C is pointed. 

Let us relabel vectors in F by /i, . . . , /m, with the first p — Si + ■ ■ ■ + Sk 
vectors belonging to F+, and denote by cyj{v) the coordinates of the vector 
V in its expansion along F. 

By (9), for v ^ Af' ■ ■ ■ A^^-f,, j = 1, . . . ,p, we have 

aj{v) > and ai{v) — for all i ^ j. 
Consequently, for such v 

p m 

J]a»>5]K»|. (11) 

j=i p+i 

Inequahty (11) obviously holds for i> e f + F^ or / + Fc, since then the 
first p coordinates (yj{v) and exactly one of the other m — p coordinates are 
equal to one, while the remaining ones are all zeros. The construction of the 
subspace L_|_ (for which F_|_ is a basis) guarantees that inequality (11) persists 
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for vectors v being images of / + F_ under arbitrary products A^^ ■ ■ ■ ^4^" . 
Indeed, the left hand side of (11) is 

EA-.-.A-", (12) 

1=1 

while the right hand side is just one summand of the form 

|A----A-"|, (13) 

with j between p + 1 and m. Since all summands in (12) are non-negative, 
and at least one of them is bigger than or equal to (13) — this is where 
the definition of L+ is being used, — inequality (11) will hold for such v. 
Moreover, for images of / + Fc under A^^ ■ ■ ■ we have, due to (10): 

P ^ m 

E«j(^) ^ g^'"^'^^'"' ^Ofor j = l,...,p. (14) 

j=i p+i 

Since inequalities (11) and (14) persist under taking linear combinations with 
non- negative coefficients and passing to limits, we see that (11) holds in fact 
for all f G /C. On the other hand, if (11) holds after switching from v to —v, 
then aj{v) — for all j = 1, . . . , m, so that v — 0. □ 



6. Families of matrices with common dominant eigenvector 

Theorem 6 gives a full treatment of families of 2 x 2 matrices sharing 
a dominant eigenvector. In higher dimensions, however, we have to impose 
additional restrictions. 

Theorem 13. Let A be a set of m x m Vandergraft matrices that share a 
common dominant eigenvector x and satisfy at least one of the following two 
conditions: 

(1) The matrices in A are simultaneously similar, with a real similarity 
matrix, to normal matrices; 

(2) A is finite, the matrices in A commute and for every A E A, p{A) is 
a semisimple eigenvalue, i.e., a simple root of the minimal polynomial, 
of A. 
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Then the matrices in A have a common invariant proper cone K, with the 
additional property that x belongs to the interior of K,. 

For the proof of Theorem 13 we need two lemmas. 

Lemma 14. Let Ai, . . . ,Aq be commuting mxm real matrices. Assume that 
there exists Xq real with the following properties: 

(1) there exists a nonzero x such that AjX = Xqx for j = 1,2, . . . ,q. 

(2) Ao is a semisimple eigenvalue of Aj, for j — 1,2, . . . ,q. 

Then there exists an invertible real matrix S such that S~^AjS have the form 



S-^A,S 



Ao 
B, 



j = 1,2, 



where Bi, . . . ,Bq are (m — 1) x (m — 1) matrices. 

Proof. Induction on q. For q — 1, the result is clear. Assume Lemma 14 
has been proved for g — 1 matrices. Applying a simultaneous similarity to 
Ai, . . . , Aq, we may assume that 



Ai 



Aoip 
Ai 



where Aq is not an eigenvalue ol A^. Since Ai, . . . ,Aq commute we have 



A, 



Bj 

a- 



j = 2,3,...,q. 



Here the matrices B2,...,Bq are p x p. Clearly, the vector x (which exists 



by (1)) has the form x 



y 





where y has p components. Then 



BjU = Xqij. One verifies that Aq is a semisimple eigenvalue of each Bj. By 
the induction hypothesis, there exists an invertible real T such that 



T-^BjT 



Ao 
B, 



j ^2,3,...,q. 



Now take S 



T 
/ 



to satisfy the lemma. 



□ 
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Lemma 15. Let Ai, . . . ,Aq be commuting m xm complex matrices with the 
following properties: 

(1) p{Aj) <1 for 1,2,..., q; 

(2) every eigenvalue (if exists) on the unit circle of every Aj is semisimple. 
Then there exists a positive definite matrix V such that 

V-A*VAj>0, for j = l,2,...,5. (15) 

(A> B means that A — B is positive semidefinite) . 

Moreover, if all Aj 's are real, then V can be also chosen real. 

Proof. It is enough to prove the complex case only. Indeed, suppose all A^'s 
are real and we have proved that there exists a (generally, complex) positive 
definite V such that (15) holds. Then by taking complex conjugates in (15) 
we obtain 

V-AjVAj>0, j = l,2,...,j = l,2,...,q. (16) 

Adding (15) and (16) we see that U - AjUAj > 0, where U -.^ V + V is 
positive definite and real. 

We now prove the complex case. If p{Aj) < 1 for all j, let 

where the sum is taken over all g-tuples {zi, . . . , Zq), zj e Z+. It is easy to 
see (using p{Aj) < 1) that the series in (17) converges absolutely. Clearly 
V>I and 

V - A*VAj 

= Y^iAir ■ ■ ■ {A*_,p-^ {A*^^P^^ ■ ■ ■ {AlpA^q^ ■ ■ ■ A^A^:; ■ ■ ■ 
> 0, 

where the sum is taken over all (g — l)-tuples {zi, . . . , Zj^i, Zj+i, . . . ,Zq) G 

So suppose that p(Aj) = 1 for some j, say p{Ai) = 1. Note that the hy- 
potheses and the conclusions of Lemma 15 are invariant under simultaneous 
similarity of Ai, . . . , Ag-. 

Aj^S-'AjS, j^l,2,...,q. 
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where S is any invertible m x m matrix. Then, considering each root sub- 
space of Ai separately, and taking advantage of the commutativity property 
AjAk = AkAj for j,k = 1,2, ... ,q, we reduce the proof to the case Ai = XI, 
|A| = 1. Then obviously V — AlVAi = 0, and it suffices to prove (15) for 
A2,. . . ,Aq. This follows by induction on q, the case q — 1 being easy. □ 

We now proceed with the proof of Theorem 13. 

Proof. Assume first that (1) holds. We may assume that A consists of normal 
matrices and that = 1 (the norm is Euclidean). Let M be the orthogonal 
complement to Span{x}. We claim that: 



is a common invariant proper cone for all A & A. 

Clearly, /C is a proper cone; therefore we only have to show that it is 
invariant with respect to the matrices. Let A E A, and let x,U2, ■ ■ ■ ,Um be 
an orthonormal set with the following properties: 

Au2k = oikU2k + PkU2k+i, Au2k+i = -PkU2k + oikU2k+i, for k ^ 1,2, . . . ,£, 

Au, = XsUs for s = 2£ + 2, 2£ + 3, . . . , m, 
where ak, (3k, are real numbers such that 



here f is a certain nonnegative integer. (The existence of such U2, . . . ,Um 
follows from the canonical form of real normal matrices, see, e.g., [13].) Ob- 
viously U2, . . . , Um form an orthonormal basis in M. Take 



;C := {cix + y : ci e M, y e M, c^>\\y\\} 




y ^CiX + C2U2 + ... + CmUm £ JC, 



thus ci > \/ + ■ ■ ■ + c^. Then we have 



Ay = p{A)cix + C2{aiU2 + P1U3) + C3{-PiU2 + aiUa) H 



+C2£{aeU2e + PeU2e+i) + C2e+i{-PeU2e + aeU2e+i) + 

X2e+2C2e+2U2e+2 + • • • + XmCmUm ■= p{A)ciX + W. 



(18) 
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Notice that for A; = 1,2, . . . ,£ we have 



= {C2k0ik - C2k+lPk)u2k + {c2kPk + C2k+lOik)u2k+l 

and 

{C2k0ik - C2k+l/3kf + {c2kPk + C2k+lOlkf _ ("^ + l3'^){4k + 4k±2) / 2 , 2 

S C2k+C2k+l- 



p{Ay 

Thus, 

\\w/p(A)\\^ <cl+---+cl,^,+ 



)i2 

^^2^+2 2 

Cn 



, + ••• + 



p(Ay 



cl<cl+---+cl, (19) 



and it follows from (18) and (19) that Ay e /C. 

Assume now that (2) of Theorem 13 holds. Let A — {Ai, . . . , Ag}. We 
may assume that the spectral radius of each Aj is positive (if some Aj is 
nilpotent, the hypotheses of Theorem 13 (assuming (2)) imply that it is 
actually equal to the zero matrix, and can be ignored). Scaling the Aj^s we 
may further assume that p{Aj) = 1, j = 1,2, . . . ,q. By Lemma 14 we may 
assume that 

" 1 



A, 







where Bi, . . . , Bq are (m — 1) x (m — 1) matrices. By Theorem 1, the hy- 
potheses (1) and (2) of Lemma 15 are satisfied for Bi, . . . ,Bq. Thus, there 
exists a real positive definite matrix V such that 



V-BjVB.>0, j^l,2,...,q. 



(20) 



Then 





X 







: X > 0, ye 



nm— 1 



is such that y^Vy < x"^ 



is a common invariant cone for ^41, . . . , Ag. Indeed, if 



A, 



X 

y 



e /C, then 



X 




X 


. y . 
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and 

{B^yfVB^y < by (20) < y'^Vy < x\ 

and so 



Clearly, KL is topologically closed, is closed under multiplication by nonneg- 
ative real numbers, is solid and pointed, because of the positive definiteness 
of V. It remains to prove that /C is convex. Thus, let Xi,X2 > and 
yi, y2 e R"""^ be such that 

ylVyk<xl for A: = 1,2. (21) 

Then for a number a between and 1, we have: 

{ayi + {l-a)y2fV{ayi + {l-a)y2) < a'xl + {l-afxl+2a{l-a){ylVy2) < 

< a^x\ + (1 — a^x^ + 2q;(1 — ajXiX'z — {axi + (1 — Q;)a;2)^ 

(Cauchy-Schwartz inequality and (21) are used in the last step of the deriva- 
tion), and the convexity of /C is proved. □ 



7. Examples 



In this section we collect examples that illuminate concepts and results 
presented. We use the notation 





" 1 " 




" " 


ei = 





e2 = 


1 



Example 1. Two 2x2 matrices A and B with negative determinants such 
that all words in A and B are Vandergraft matrices though there is no {A, B)- 
invariant proper cone. 



Take 



1 p 
-1 



B 



1 q 
-1 



pj^q. 



All words in A and B are Vandergraft matrices, with ui = ei as a domi- 
nant eigenvector. So, Theorem 6 applies, and according to case (iii) in "Only 
if" part of its proof {A, i?)-invariant proper cones do not exist. 

Example 1 shows that Theorem 7.6 in [11] is apparently misstated. 
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Example 2. A triple of matrices T := {A, B, C}, A,B,C e Ry^^ with the 
following properties: 

(a) det M > for all M e T; 

(b) A,B,C are normal matrices (in particular, diagonalizable); 

(c) there is no T-invariant proper cone; 

(d) each pair of matrices in T has a common invariant proper cone; 

(e) no two matrices in T have a common eigenvector. 

The example shows that sharing a common dominant eigenvector is es- 
sential in Corollary 3 and Theorem 13, and also that the part of Theorem 9 
pertinent to the case when there are no simultaneously diagonalizable pairs 
of matrices in A is sharp. 

Instead of describing the matrices directly, we will list two linearly in- 
dependent eigenvectors and associated eigenvalues for each matrix. For the 
eigenvalues simply pick Ai(M) > \2{M) > for each matrix M E T. As for 
the eigenvectors of a matrix Af , denoting the dominant and non-dominant 
ones by Ui{M) and U2{M) respectively, let 



and 



ui{A) = ei, ui{B) 



U2{A) = 62, U2{B) 



U2{C) - 



1 

-2 
1 



Each of the pairs {A, B), (A, C) and {B, C) then satisfies conditions of Corol- 
lary 4, and therefore has a common invariant proper cone (more specifically, 
Cone{ui{A),ui{B)} is (A, 5) -invariant, Cone{— Mi(C)} is 
{B, C)-invariant, and Cone{ui{A) , ui{C)} is {A, C)-invariant). On the other 
hand, the separation condition (ii) of Corollary 4 does not hold for the triple 
(^4, B, C), so that there is no (^4, B, C)-invariant proper cone. 

Example 3. A quadruple of matrices A, B,C,D & R^^^ with distinct positive 
eigenvalues such that each triple of them has a common invariant proper cone 
while there is no {A, B, C, D)-invariant proper cone. 
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In accordance with Theorem 9, this quadruple consists of two pairs of 

commuting matrices. 

As in Example 2, the eigenvalues of the matrices can be chosen arbitrarily, 
as long as they are positive and distinct. Following the eigenvector notation 
from the same Example, let 

ui{A) = U2{B) = ei, U2{A) = ui{B) = 62, 

ui{C) = U2{D) = ei + 62, U2{C) = ui{D) = ei - 62- 

The vectors 61,62,61 + 62,61 — 62 are simultaneously dominant and non- 
dominant for the quadruple {A, B,C, D), and cannot be separated in the 
sense of condition (iii) of Theorem 7. Consequently, there is no {A, B, C, D)- 
invariant proper cone. On the other hand, from Corollary 4 it follows (and 
can also be checked directly, based on Theorem 2) that Conc{ei,e2} is 
{A, B, C)-invariant, Cone{6i, —62} is {A, B, -invariant, Conc{6i + 62, 61 — 
62} is (^4, C, £)) -invariant, and Cone{6i -|- 62, 62 — 61} is (S, C, £))-invariant. 

Exampl6 4. The set S = {A, B} which satisfies all the hypotheses of The- 
orem 13 (with (2) holding) except that p{A), p{B) are not semisimple eigen- 
values of A, B, respectively, and there is no {A, B) -invariant proper cone. 



Take 



A 



1 1 
1 



B 



1 -1 
1 



Clearly, both matrices arc Vandergraft, non-diagonalizable, sharing the dom- 
inant eigenline but having different orientation. By Theorem 7, there is no 
common invariant proper cone. 

Example 5. Two diagonal matrices Ai and A2 without a common invariant 
proper cone such that all words in Ai and A2 are Vandergraft matrices: 



A, 



1 
-1 









■1 



A, 









-1 
1 



It is easy to check that all words in Ai and A2 are Vandergraft matrices. 
However, condition (9) of Theorem 12 fails, so that there is no (^41,742)- 
invariant proper cone. 
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Example 6. Countable set 0/ 2 x 2 Vandergraft matrices such that every 
finite number of them has a common invariant proper cone, but the whole set 
does not- 
Using Theorem 2, it is easy to see that any set of the form 



where the sequence {|?m|}m=i tends to infinity and < r < 1 is fixed, fits 
the bill 

Remark 3. From standard compactness considerations it follows that if A 
is an infinite family in M^^" any finite subfamily of which has a common 
invariant proper cone, then there exists a non-trivial (that is, different from 
{0}) A-invariant closed convex pointed cone. However, it may not be solid, 
and therefore is not necessarily proper. 

This is exactly what is happening in Example 6. 
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